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Method, apparatus and software for analysing perfusion images 



The invention relates to a method for analysing perfusion images, in particular 
MR perfusion images, of a human or animal organ including the steps of 

(a) defining at least one contour of the organ, and 

(b) establishing at least one perfusion parameter of a region of interest of said 
5 organ within a boundary defined by the at least one contour. 

From the article fl MR perfusion imaging: correlation with PET and 
quantitative angiography", published in Magnetic Resonance Materials in Physics, Biology 
and Medicine 1 1 (2000) 71-72 by J. Schwitter en G.K. von Schulthess, it is known to apply 
quantitative tools for obtaining pixelwise slope maps of registered MR perfusion data sets in 
10 a method for detecting stenosed coronary arteries from the heart. Analysis of upslope 

(meaning perfusion rates), rather than other parameters is said to provide a very sensitive and 
specific measure of myocardial ischemia. Measuring the upslope provides a semiquantitative 
measure of absolute perfusion, even with patients having triple vessel disease. 

The invention aims at improving the method mentioned in the introductory 

15 paragraph. 

To this end, the method of the invention is characterised in that steps (a) and 
(b) are repeated in a series of iterative steps wherein for each subsequent iterative step the 
definition of the at least one contour in step (a) is varied, and the series of iterative steps is 
terminated after reaching an optimal value for the at least one perfusion parameter in step (b). 

20 By this measure, a very accurate determination of an organ's boundaries is 

possible allowing for a consequently also very accurate determination of the perfusion 
parameter or parameters that are to be established. 

The method of the invention is particularly useful when the organ is a heart 
and the region of interest is the heart's myocardium or a segment thereof. 

25 It is remarked, however, that the gist of the invention and the scope of 

protection afforded hereby is not limited to analysing perfusion images with respect to the 
heart's myocardium or one or more segments thereof The invention may equally well be 
applied in respect of other organs, for instance the brain, The above and the further 
explanation below with respect to the analysis of perfusion images of the heart, and in 
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particular the heart's myocardium, serves therefore primarily to elucidate the method without 
being restricted to that single application. 

preferably, in step (a) the inner contour and/or the outer contour of the heart's 
myocardium is defined. Although one could resttict the optimisation according to the method 
5 of the invention to variation of the myocardium's inner contour, particularly good results are 
attainable if both the inner contour and the outer contour of the heart's myocardium are 
varied. This variation may occur simultaneously, but it is preferred to execute the variations 
subsequently. 

Suitably the perfusion parameter to be established is selected from the group 
10 comprising mean upslope, time to peak intensity, time of contrast arrival, tune to half 
maximum intensity, accumulated inflow or combinations thereof. It is found particularly 
advantageous however that in step (b) the perfusion rate (or upslope), and/or the time at 
which the maximum perfusion rate occurs is established. These parameters provide effective 
indication of the perfusion of the organ being analysed. 
15 in caS e the myocardium's inner contour is varied, the method of the invention 

is particularly well applied by executing the series of iterative steps and ter mina t in g these 
series after reaching an essentially constant value for the perfusion rate and/or said time at 
which the maximum perfusion rate occurs, as compared to the perfusion rate's value and/or 
time in a previous iterative step. 
20 The invention is further embodied in a software program for a computer of an 

apparatus that is implemented to execute the method for analysing perfusion images, as 
explained above. 

The invention is further embodied in an apparatus that is apparently intended 
to execute said method. 

25 The invention shall now be further elucidated with reference to the following 

non-limiting exemplary embodiment of the invention with reference to the drawing. 
In the drawing it is shown 
in Fig. I : a short axis slice through a heart; 

in Fig. 2: intensity time curves of a bolus passage through sections of the 

30 heart; 

in Fig. 3: perfusion parameters varying depending on radial displacement of an 

inner myocardial boundary; and 

in Fig. 4: variations of perfusion parameters depending on displacements of a 

myocardial outer boundary, 
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In the example, the method for analysing the perfusion images relates to the 
myocardium 1 of the left ventricle 2 of a heart (see Fig. 1), Reduced blood perfusion of the 
myocardium 1 of the left ventricle 2 is a direct result of cardiovascular diseases. For 
measuring the perfusion of the myocardium 1, a contrast fluid is applied whilst the patients 
5 heart is monitored with MR or other imaging techniques that are known per se. The images 
are obtained by making scans of slices 3 through the myocardium 1 during a period of 20 to 
40 seconds. In an image sequence, it can be monitored that the contrast fluid first enters the 
right ventricle 4, then the left ventricle 2, and finally the myocardium 1. 

The perfusion of these parts of the heart is imaged and the imaging is 
10 supported by the contrast agent. The intensity of the measurement as recorded in time, is 
shown in Fig. 2. The X- axis of Fig. 2 shows time as an independent variable, whilst the Y- 
axis shows the intensity of the measured perfusion* The Figure clearly shows that the contrast 
agent first enters the right ventricle. This is shown by graph A. Subsequently, the contrast 
agent arrives at the left ventricle as shown by graph B, Finally the contrast agent arrives at 
15 the myocardium as shown by graph C. The intensity time curve C pertaining to the 
myocardium is used for analysing the perfusion of the myocardium. 

Since "die signal in the myocardium 1 is very noisy, generally the myocardium 
1 is divided into segments and the measurements are averaged over these segments. These 
segments preferably coincide with areas of the myocardium 1 which are supplied with blood 
20 from a certain coronary artery. In this way, if a reduced perfusion is observed in a myocardial 
segment, it can be traced back to the supplying artery. 

Since the signal strengths in the left 2 and right 4 ventricles are much higher 
than in the myocardium 1 , extreme care should be taken not to include any part of the left and 
right ventricles into the myocardial segments. This is not as easy as it seems, since the 
25 boundaries of the left and ventricle blood volumes can be very irregular due to the presence 
of papilar muscles and trabeculae. 

Also the boundaries of the myocardium 1 are often not clearly visible in a 
single image of the sequence or even on a maximum intensity projection through time. 
Furthermore, the images generally have a rather low resolution (typically 128*128 pixels) 
30 and only a few pixels are available for averaging within a segment 

The parameters that are of interest in analysing the perfusion are on the one 
hand the perfusion rate or upslope as indicated in Fig, 2 with line D in relation tp graph C, 
and on the other hand the time t u at which the maximum perfusion rate occurs. Fig. 2 relates 
to a specific situation with a fixed inner boundary and outer boundary of the myocardium 1 . 
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Figs. 3 and 4 show the results with a varying diameter of inner contour and 
outer contour of myocardium 1 respectively, and the resulting value of the maximum upslope 
or maximum perfusion rate, and the time at which this maximum perfusion rate occurs (see 
Figs. 3 and 4 left and right respectively). 
5 The X-axis of Figs. 3 and 4 shows as the independent variable a variation of 

the concerning contour of the myocardium (Fig. 3: inner contour; Fig. 4: outer contour). 

From Fig. 3 is follows that when the myocardium's inner contour is varied, 
eventually the above-mentioned perfusion parameters reach an essentially constant value, 
whilst further variation of the inner contour, i.e, a diminidiing diameter beyond the point at 
1 0 which the constant value is reached, does no longer effect either the maximum perfusion rate 
or the time at which such maximum perfusion rate occurs. The point beyond which further 
variation of the inner contour does no longer effect the perfusion parameters is then 
estahlished to represent the measured inner contour that most accurately corresponds to the 
true boundary of the myocardium. 
1 5 Likewise* when varying the outer contour of the myocardium, an optimal 

value of both the perfusion rate and the time at which the maximum perfusion rate occurs, 
can be established (see Fig. 4). The outer contour value corresponding to these optimal values 
is then taken to accurately represent the actual outer boundary of the myocardium. 

When one or more optimal values of the said perfusion parameters are 
20 established, further iteration by repeatedly defining a contour of the organ and subsequently 
establishing at least one perfusion parameter of the analysed region of interest of the organ 
can be terminated. In this way, an accurate method for detecting the myocardial contours is 
defined that ensures that no part of the left ventricle 2 or right ventricle 4 is included in the 
perfusion images of the myocardium 1. 
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CLAIMS: 



10 



1 . Method for analysing perfusion images, in particular MR perfusion images, of 
a human or animal organ including the steps of 

(a) defining at least one contour of the organ, and 

(b) establishing at least one perfusion parameter of a region of interest of said 
organ within a boundary defined by the at least one contour, 

characterized in that steps (a) and (b) are repeated in a series of iterative steps wherein for 
each subsequent iterative step the definition of the at least one contour in step (a) is varied, 
and the series of iterative steps is terminated after reaching an optimal value for the at least 
one perfusion parameter in step (b). 

2. Method according to claim 1, characterized in that the organ is a heart and the 
region of interest is the heart's myocardium or a segment thereof. 



3. Method according to claim 2, characterized in that in step (a) the inner contour 
1 5 and/or the outer contour of the heart's myocardium is defined. 

4. Method according to any one of claims 1-3, characterized in that in step (b) the 
perfusion rate or upslope and/or the time at which the maximum perfusion rate occurs is 
established. 

20 

5. Method according to claim 4, characterized in that in step (b) the 
myocardium's inner contour is varied and that the series of iterative steps is terminated after 
reaching an essentially constant value for the perfusion rate and/or said time at which the 
maximum perfusion rate occurs* as compared to the perfusion rate's value and/or time in a 

25 previous iterative step. 

6. Software program for a computer of an apparatus implemented to execute the 
method according to any one of claims 1-5. 
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7, Apparatus comprising mean9 to execute the method according to any one of 

claims 1-5. 
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ABSTRACT: 



The invention relates to a method for analysing perfusion images, in particular 
MR perfusion images, of a human or animal organ including the steps of 

(a) defining at least one contour of the organ, and 

(b) establishing at least one perfusion parameter of a region of interest of said 
5 organ within a boundary defined by the at least one contour, whereby steps (a) and (b) axe 

repeated in a series of iterative steps wherein for each subsequent iterative step the definition 
of the at least one contour in step (a) is varied, and the series of iterative steps is terminated 
after reaching an optimal value for the at least one perfusion parameter in step (b). 
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